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Abstract 

Broadening of the Landau levels in graphene and the associated quantum Hall plateau-to-plateau transition are investigated numer- 
ically. For correlated bond disorder, the graphene-specific n - Landau level of the Dirac fermions becomes anomalously sharp 
accompanied by the Hall transition exhibiting a fixed-point-like criticality. Similarly anomalous behavior for the « = Landau level 
is also shown to occur in correlated random magnetic fields, which suggests that the anomaly is generic to disorders that preserve 
the chiral symmetry. 
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1. Introduction 



After the seminal discovery of the characteristic quantum 



Hall effect (QHE) in graphene |[1 , 
on the effect of disorder Jsl 0, 



0, attention has been focused 
Randomness, which is 



important in the QHE in ordinary 2DEG, is particularly cru- 
cial since it affects such key factors in the lattice (as opposed 
to Dirac field) models of graphene such as the chiral (A-B sub- 
lattice) symmetry 1 8], and the scattering between the valleys (K 
3nd K'). Note that the valley degeneracy has a topological ori- 
gin as the doubling of the massless Dirac fermions on the hon- 
eycomb lattice. It is then natural that the nature of disorder 
crucially affects the electronic structure. A potential disorder 
(including random site energies) destroys the chiral symmetry, 
whereas a disorder in bonds respects the symmetry, while the 
spatial correlation in the disorder controls the inter-valley scat- 
tering. In order to explore these, here we adopt the honeycomb 
lattice rather than the effective Dirac model to investigate how 
the spatial correlation of disorder affects the Landau level struc- 
ture, especially the stability of the n = Landau level which is 
essential to the characteristic QHE in graphene. 

In the case of a potential disorder with degraded chiral sym- 
metry, the quantum Hall transition at the n = Landau level 
has been shown to be robust in an effective Dirac model |0] . In 
such a model, however, the criticality of the transition at n = 
is described by the ordinary quantum Hall transition, with noth- 
ing special about the criticality at the n - transition. In the 
previous paper [9], we have shown that for the bond disorder, 
which preserves the chiral symmetry, the criticality at « = 
transition is anomalously sensitive to the spatial correlation of 
disorder. As soon as the correlation length of the bond disor- 
der exceeds few lattice constants of the honeycomb lattice, the 
n - Landau level becomes anomalously sharp and the asso- 
ciated Hall transition exhibits the anomalous criticality, which 



corresponds to a fixed point of the Dirac fermion with chiral 
symmetry lilQI . In the present paper, we present another exam- 
ple to reinforce our arguments that the anomalous Hall transi- 
tion is a general property of the chiral-symmetric random sys- 
tems. Namely, we consider systems with random phases in the 
transfer integrals, which corresponding to the case of random 
magnetic fields piercing the hexagons in honeycomb lattice. 

2. Models 

The tight-binding model for the honeycomb lattice is de- 
scribed by the Hamiltonian // - tije^^''^^'Jclcj-i-h.c., where 
the transfer integral tij is real while the Peierls phases Otj satisfy 
the requirement that the sum of the phases around a hexagon is 
equal to the magnetic flux piercing a hexagon in units of the 
flux quantum tpo - h/e. 



Random bonds are introduced as f. 



t H- dtij, where t and 
6tjj are the uniform and the random components, respectively. 
The random components 6tij are assumed to have a Gaussian 
. g-(S/-/2a-2 1 ^2ncrj, with a variance cr, and a 



distribution, P(6t) ■ 
spatial correlation. 



{dtijdhi) = {6t^)exp(-\rij-rkif-/4-r]j), 

with a correlation length 77,, where denotes the position of 
the bond and () the ensemble average. 

On the other hand, random magnetic fields are introduced as 
0(r) - <p -i- 6tp{r). This type of disorder, being another disorder 
in bonds, also preserves the chiral symmetry. Here ^ represents 
the uniform part of the magnetic field, while 6(p(r) the random 
magnetic fluxes for hexagons each located at position r. We 
assume that the random component 6(p obeys a Gaussian distri- 
bution with a variance cr^ and a spatial correlation 



{6(l){ri)6<l){rj)) 



-} exp(- 



7H) 
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with a correlation length T/^fTT]. In the following, all lengths 
are measured in units of the bond length a of the honeycomb 
lattice. 

Let us recapitulate how these types of disorder preserve the 
chiral symmetry (even for each realization of disorder). There 
exists a local unitary operator y that anti-commutes with the 
Hamiltonian, {y,H] - 0, which defines the chiral symmetry, 
and = 1. In the case of the honeycomb lattice, the lattice 
sites can be decomposed into A and B sublattices . The oper- 
ator y can then be given as y = exp(i7r Y^ieA c'^c,), where YtieA 
denotes the summation over the A sublattice sites. With this 
choice of y we can readily verify that the fermion operator c, 



is transformed as yciy 



-Cj for / e A and ycjy 



for i e B. Due to this chiral symmetry, the energy levels al- 
ways appear in pairs {E, -E} even in the presence of disor- 
der For instance, if we have an energy eigenstate if/E satisfying 
HtJ/E - Eif/E, the state yij/E is an eigenstate with energy -E 
since Hytf/E = -yHiffE = -EyiffE- A special situation arises for 
zero-energy states, for which the eigenstates 4'e=o and yi/'£=o 
are degenerated. The eigenstates can then be made simulta- 
neous eigenstates i/r-t = (1 ± y)4'E=Q of the operator y with 
yili± - +ifr±). The eigenstate has non-zero amplitudes 

only on the A{B) sublattice. The fact that the zero-energy states 
can be an eigenstate of y implies that the zero-energy states has 
an extra symmetry, so that the criticality at zero energy can be 
especially sensitive to the presence or otherwise of the chiral 
symmetry \12, T^. Obviously, the potential disorder breaks 
this symmetry. 

The density of states (p,) = - 2,- ImGuiE + \e)INn is evalu- 
ated by the Green function iH] GuiE + ie) = {i\{E-H +ie)-^\i), 
where stands for the total number of sites and e an infinites- 
imal imaginary part of energy for evaluating the Green func- 
tion numerically. We have carried out calculations for values 
of e/f reduced from 0.01 down to 6.25 x 10""* to confirm that 
the anomaly in the density of states at the n - Q Landau level 
discussed below is not affected by the value of e. The system 
size considered is a x Ly rectangular system with periodic 
boundaries in y direction, where x axis is assumed to be paral- 
lel to the zigzag direction of the honeycomb lattice. In actual 
calculations, the Landau gauge for the corresponding bricklayer 
lattice lll5r is adopted. 



3. Numerical Results — Random Bonds 

First, we consider the random bonds in a uniform magnetic 
field ((7-0 = 0). In the previous paper we have clearly 
demonstrated for such a model that when the correlation length 
T], exceeds few bond lengths, the n = Landau level becomes 
anomalously sharp and the associated quantum Hall transition 
shows an almost exact fixed-point-like criticality even in a finite 
system (Fig. [TJ. This sharply contrasts with the result for the 
uncorrected case = , for which the n - Q Landau level is 
broadened in the same way as n Landau levels. The anoma- 
lous sensitivity to the correlation of disorder occurs only for the 
n -Q Landau level. We have examined the density of states for 
other values of o", and 0, and confirmed that the anomaly at the 
n = Landau level is commonly observed for rjja > 1. Here 




E/t 



Figure 1: The averaged Chem number nc (the Hall conductivity {cr yy / (e-^ / h))) 
as a function of the Fermi energy (E/t). The Chern number for each realization 
of disorder is evaluated by using the lattice gauge technique(l,[ll[li|Tl. 
The results for the correlation length r],/a = 0, 0.5, 1.0, and 1.5 are shown for 
plateaus nc = <Tyy,/(e-/h) = -3, -1, 1, 3 with the disorder strength a-,/t = 0.12, 
the magnetic field (p/cpo = 1/50, and an average over 300 samples for the system 
size LJi a/3£!/2) = Ly/(3a/2) = 20. 



we show the density of states as a function of the Fermi energy 
and the correlation length rj, for cr,/r = 0.058 and 0/(^o = 1/50 
in Fig. |2] It is again clearly seen, as has been seen in the case 
of cTi/t - 0.12 CT, that the n - Landau level becomes anoma- 
lously sharp compared to other Landau levels as soon as the 
correlation length rj, is greater than the lattice constant. 

In actual graphene samples, the scale of ripples is estimated 
to be of order of 10 nanometers 1, 19. ,201 . which is much greater 
than the lattice constant a ~ 1.42 A ||7[], and the magnitude 
of disorder is likely to be much smaller Our result there- 
fore clearly indicates that the bond disorder induced by ripples 
should not broaden the « = Landau level. 



4. Numerical Results — Random Magnetic Fields 

Let us next show the results for another model, namely the 
random magnetic field model with cr^ 0, where no ran- 
domness is assumed for the amplitude of the transfer energies 
(cr, - 0). For this model, we consider both the small random 
field case, cr^ < (p (case 1), and the large random field case, 
cr^> <p (case 2). For cr^ < (p the distinct Landau level structure 
is present (Fig. [3]l. The n = Landau level is again anoma- 
lously sensitive to the spatial correlation of random magnetic 
fields: When the correlation length tj^ exceeds few lattice con- 
stants, the width of the n - Q Landau level becomes a sharp, 
delta-function-like peak. Indeed, the shape of the n - Landau 
level for rj^la > 2.5 is almost exactly given by the Lorentzian 
distribution (<p/(po){£/7:)/(E^ + e^), indicative of a zero intrinsic 
width. This again contrasts with the « Landau level. 

For (Tip > ^ (case 2), on the other hand, the Landau level 
structure is mostly washed out for n + Q Landau levels, which 
is not surprising since the fluctuation of the magnetic field is 
larger than its mean value. Surprisingly, n - Q Landau level 
exhibits an anomalous behavior, where the delta-function-like 
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Random Bonds 




Figure 2: The density of states for various values of the correlation length Tj,la 
of the random bonds. The system-size is LJi^alT) = 5000, Ly/(3a/2) = 
100. The imaginary part of energy e/t = 6.25 X lO"'*. 

behavior arises when the correlation length of the random mag- 
netic field exceed few lattice constants, even though the field 
fluctuation width is much greater than the average field (FiglU. 

5. Discussions and Conclusions 

In order to check that the above anomalies are a conse- 
quence of the preserved chiral symmetry, we have also eval- 
uated the density of states for the case of the spatially corre- 
lated potential disorder, where the Hamiltonian is given hy H = 
'Zj(ij)(te^^''^'^c'. Cj+h.c.)+Yji sic\ Ci, where random site energies s,- 
are assumed to obey a Gaussian distribution with a variance cr^ 
and a spatial correlation (e/fij) = cr^ exp(-(r/ - rjf' jAif;.). The 
result for cr^/f - 0.29 and (pl(pQ - 1/41 in Fig. |5]confirms that 
no anomaly exists at the n - Landau level. As we increase 
the correlation length rjs keeping the magnitude of disorder ctj 
fixed, all the Landau levels are broadened [21i] . 

To summarize, we have investigated the anomalous behav- 
ior of the n = Landau level of the QHE system on a honey- 
comb lattice. Two types of bond disorder that preserve the chi- 
ral symmetry have been considered, namely the randomness in 
the magnitude of the transfer integrals (random bonds) and that 
in their phases (random magnetic fields). It is clearly demon- 
strated that in both cases the n = Landau level is anoma- 
lously sensitive to the spatial correlation of disorder, where an 
anomaly at the n - Q Landau level appears as soon as the corre- 
lation length exceeds few bond lengths. This indicates that the 
absence of the mixing between two valleys are essential to the 
anomaly at the n - Q Landau level. We have also confirmed that 
the anomaly does not exist for the correlated potential disorder, 
where the chiral symmetry is broken. The results suggest that 
the anomaly at the n - Q Landau level is generic to the graphene 
system with long-range disorders that preserve the chiral sym- 
metry. This implies that the bond disorder induced by ripples 



Random B < Uniform B (Case 1) 




Figure 3: The density of states is plotted for various values of the con'elation 
length rjif, /a when the fluctuation of the magnetic field is smaller than its mean 
(o-0 < (p). The parameters are a-^/cpo = 0.0058 < (p/(po = 1/41 - 0.024 for 
a system size Lv/( V3«/2) = 5000, Ly/(3a/2) = 82 and the imaginary part of 
energy e/f = 6.25 X lO""*. 



Random B > Uniform B (Case 2) 




Figure 4: As above when the fluctuation of the magnetic field is greater than its 
mean (cr^ > </>). The parameters are cr^/ipQ = 0.058 > 0/0o = 1/41 - 0.024 
{(/> < o-^) for a system size L,-/(V3a/2) = 5000, L,7(3fl/2) = 82 and the 
imaginary part of energy e/f = 6.25 X lO"*. 
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Figure 5: The density of states in a random potential for various values of the 
coiTelation length r]sla. The parameters are (Tj/( = 0.29, <l>l(j>a = 1/41 and 
0-, = O-0 = for a system size L^K^all) = 5000, Ly/(3a/2) = 82 and the 
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in graphene should not broaden the n - Landau level; con- 
versely, if a broadening is observed in experiments, that must 
be caused by other origins, such as the potential disorder from 
charged impurities in substrates. Experimentally, a narrower 
« = Landau level of grephene has been reported by measur- 
ing the activation energy gaps ll22ll . which may be related to the 
present analysis. More elaborate analysis on the nature of the 
disorder is an interesting future problem. 
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